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Abstract 



Recent interest in the study of color superconductivity has focused on the 
regime of high baryon density where perturbative QCD may be employed. 
Based on the dominant one-gluon-exchange interaction, both the transition 
temperature and zero temperature gap have been determined to leading order 
in the coupling, g. While the leading non-BCS behavior, Tq ~ fig~ 5 e~ K ^ 9 , 
is easily obtained, the pre-exponential factor has proved more difficult to 
evaluate. Focusing on the transition temperature, we present a perturbative 
derivation of this factor, exact to leading order in g. This approach is first 
motivated by the study of a toy model and involves working to second order 
in the perturbative expansion. We compare this result to the zero tempera- 
ture gap. Additionally, we extend the analysis to the case of higher angular 
momentum for longitudinal and transverse quark pairing. 
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I. INTRODUCTION. 



Previous studies of the quark-gluon plasma at high baryon density have indicated the 
possibility of a pairing instability which may induce a novel superconducting phase of QCD 
Jl|-Q|. Based on conventional superconductivity theory, it is known that even the weakest 
attractive interaction between quarks would be sufficient to form a superconducting pair; 
for SU(3) C QCD this is provided by the 3 channel in the diquark interaction. Motivated 
by studies of QCD under unusual conditions, interest in this phenomenon has recently 
been revived @||. In fact the phase structure of matter in this regime has important 
physical consequences for the behaviour of dense astrophysical objects as well as for heavy 
ion collisions and has raised fundamental questions about the nature of QCD. 

In general, for SU(N) C QCD with Nf light flavors, the combined symmetry of the theory 
is initially SU(N) C x SU(Nf)i, x SU(Nf)R x U{1)b- However, a diquark condensate, (qq), 
would necessarily carry both color and flavor charge, and hence break the above symmetry. 
The pattern of symmetry breaking is especially interesting for the case of three flavours and 
colors, in which case the initial SU(3) C x SU(3)l x SU(3)r x U{1)b symmetry is broken to 
the global diagonal subgroup SU(3) c+ l+r P-fLO]. Since the condensate remains symmetric 
under combined color and flavor rotations, it "locks" the color and flavor degrees of freedom. 
Further interest in the color-flavor locked state has arisen from the observation that many 
properties such as quantum numbers and roughly similar masses are shared by the bosons 
produced by this symmetry breaking and hadronic matter, suggesting a connection between 
these two phases of QCD 0. These results have led to further investigation of the color- flavor 
locked phase fill and the development of an effective field theory for the Nambu-Goldstone 
bosons produced by the symmetry breaking |12| |T~T| . 

The nature of the phase transition to the quark gluon plasma and the hadronic phases 
of QCD has also been studied. Due to competition between confinement and chiral conden- 
sation on one side and the formation of a diquark condensate on the other, it is expected 
that the transition to the hadronic phase is of first order [f^pl5|-|T8|[ . Furthermore the phase 



transition to the quark gluon plasma has also been identified as first order [15,18]. Other 



physical implications that have been studied include the effect of differing densities of s 



compared to u and d quarks, [19], and the effect of a color-superconducting core upon the 



magnetic field of neutron stars, observed to be stable over very long time scales p]|-|2"5|] . 

Since the study of color superconductivity involves working at non-zero chemical poten- 
tial, it has not been amenable to numerical simulations on a lattice. Instead, a number 
of analytic methods have been developed to qualitatively probe the color superconducting 
state @-§|,||JT^,|T6|j2^-^] . All of these approaches, to some extent, have difficulties with pos- 
sible Ansatz dependence, gauge invariance, the Meissner effect as well as the non-linearity 
of the gap equation, thus making calculations of the condensate properties or the transition 
temperature somewhat difficult. The consensus of all these models, however, is that, as 
expected, the strength of the quark-quark interaction in the most attractive channel is of 
fundamental importance in determining the size of the condensate and hence the properties 
of the superconducting state |L8] . 



Quantitative results on color superconductivity may be obtained in the limit of ultra- 
high baryon density, where perturbative QCD is expected to be valid. In this regime, the 
dominant attractive interaction arises by single gluon exchange and may be incorporated in 
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a Schwinger-Dyson analysis to yield either superconducting gap equations or conditions on 
the transition temperature. Because of the long range nature of the single gluon exchange 
interaction, the coupling dependence of the gap A takes on the non-BCS form 
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and this behavior holds for the transition temperature Tq as well. In this paper we refine 
the analysis of the transition temperature and present an exact calculation of Tc to leading 
order in the coupling of the exponential (n) and pre-exponential (c) factors. This is carried 
out for both longitudinal (LL and RR) and transverse (LR) quark pairings of arbitrary 
angular momentum, J. We find that the s-wave (J = 0) channel, in which only longitudinal 
pairings of quarks may occur, is dominant. For SU(N) C QCD the general result takes the 
form, 
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where Nf is the number of flavors, fi is the chemical potential and g is the running coupling 
constant evaluated at \x. The largest factor, 



An = 256tt 4 
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was first calculated in p8|- p0|1 , initially using an elegant renormalisation group analysis. The 
pre-factors C 1 and represent angular momentum channel independent and dependent 
corrections respectively, the calculation of which are the focus of this paper. This novel 
non-BCS dependence of An on the coupling has been confirmed in H3TH331 . In addition, 
the leading order ratio of the transition temperature to the zero temperature gap (in the 
s-wave channel) was investigated in |33| and found to be identical to the BCS case, namely 
knT = —A where 7 = 0.5772 ... is the Euler constant. 

In order to exactly determine the pre-factors C l and C$ , which are of 0(1), it is most 
straightforward to approach the transition temperature from the normal phase rather than 
from within the superphase. In the normal phase there is no condensate, no Meissner effect, 
the problem is linear and there is an established framework to ensure gauge invariance. 



Following the formulation developed in [34] for the non-relativistic Fermi gas, the Dyson- 
Schwinger equation for the scattering of two quarks can be cast into a general eigenvalue 
problem. This is outlined in |3]J and more details are given in this paper. The general 



eigenvalue problem can be solved to give a set of eigenvalues, {A^}, in terms of the parameters 
of the theory. The transition to the superphase is then characterised by the locus on the 
phase diagram where the smallest eigenvalue becomes one. Physically, this condition defines 
the pole in the solution to the integral equation which represents the formation of a bound 
state of quarks. 

The pre-factor C\ is found to be of 0(1), and has the form: 
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To understand the physical origins of this factor, it can be decomposed into three parts, 
C\ = d x d[d" y each of which are entirely due to single gluon exchange. Higher order diagrams 
have been shown not to contribute to the leading order of the exponential or the pre-factor in 
The first contribution, c[, originates in the radiative corrections to the quark propagator 
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due to its self-energy [31 
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The remaining two contributions, c'[ and c'f, are related to the next-to-leading order logarith- 
mic infrared behaviour overlooked in the obtention of ( |1.3| ). The factor c'{ accounts for the 
sub-leading behaviour lost in the transition from the Matsubara sum at finite temperature 
to a continuous energy integral, 



2e^ 
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The final factor, c'(', represents the contribution of all other approximations made in the 
obtention of ( |1.3| ). Perhaps somewhat surprisingly, we find no such contributions to leading 
order in the coupling, so that in fact 
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The second pre-factor, C$, classifies the dependence of the critical temperature upon the 
angular momentum of the pairing. The dependence upon angular momentum is different for 
longitudinal and transverse condensates. In the dominant s-wave channel, for longitudinal 
pairs we find O/ =0 = 1 whereas, due to conservation of angular momentum, for transverse 
pairs we find O/ = ° = 0- Both longitudinal and transverse pairs may exist for higher angular 
momentum and the results are summarised below, 
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longitudinal pairs, J > 0, 
transverse pairs, J > 1, 



where, 



cj 



for J = 0, 

-2E n J =1 J forJ^O. 
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This indicates the suppression of the higher angular-momentum channels, so that s-wave 
pairing is dominant and p-wave pairing is already down by about three orders of magnitude. 

Before proceeding to the complete derivation of fll.2|) , we shall first illustrate the per- 
turbative formalism in section II and apply it to an exactly soluble toy model. Then in 
section III, the Dyson-Schwinger equation for the interaction between a longitudinal pair of 
quarks of zero angular momentum is studied in the Coulomb gauge. The resulting Fred- 
holm equation is cast into a general eigenvalue problem. A perturbative solution for this 
eigenvalue problem is found in section IV. Systematic calculation to second order in this 
perturbation series captures all of the double and linear logarithmic infrared behaviour of 
the eigenvalue that determines ( |1.2|) to leading order in the coupling. The generalization 
to non-zero angular momentum and transverse pairings is made in section V and in section 
VI we discuss numerical verification and the gauge invariance of our results. Concluding 
remarks are made in section VII. 
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II. A TOY MODEL: THE TWO-DIMENSIONAL CYLINDRICAL WELL. 



In this section we shall develop the perturbative formalism that will be used in the 
following sections. As shown in |j31| , the Schwinger-Dyson equation determining the tran- 
sition temperature has the form of a Fredholm integral equation and may be re-cast as 
an eigenvalue problem. We thus seek a systematic perturbative method of developing the 
eigenvalues of the Fredholm equation. To illustrate this method, we begin by considering 
the Schrodinger problem for a two-dimensional cylindrical well of depth K, 2 /a 2 and radius a: 



(where we have normalised the mass so that 2m =1). As a toy model, it is easily solved 
in terms of Bessel functions and matching conditions at the well boundary. On the other 
hand, it may also be converted into an integral equation of Fredholm type. The perturbative 
solution to the resulting Fredholm equation will serve the basis for generalization to the real 
problem of determining the QCD transition temperature. 

Note that although two parameters appear in ( |2.1| ), namely a and n, only the dimen- 
sionless coupling k is physical, as a may be scaled away by x — > ax. In analogy with 
the QCD problem, k represents the QCD coupling g, and a represents the inverse of the 
chemical potential [i. The eigenvalue A 2 plays the role of either the critical temperature 
or the mass gap of the color superconductor and in the weak coupling limit has the form 
A oc (1/a) exp(— l//t 2 ). While this behavior is BCS-like (as may be deduced from the short- 
range nature of the cylindrical well potential) some mathematical features of this toy model 
nevertheless extend to the long-range QCD case. 

In the following, we show that this leading order in k behaviour of A can be recovered 
from a perturbative expansion of the resulting Fredholm kernel, where the zeroth order op- 
erator is defined to capture the leading logarithmic behaviour. In this case we find that 
the exponential and its pre-factor can be determined to leading order in k by second order 
perturbation theory. Furthermore, the exponential is determined by the leading logarith- 
mic behaviour and the pre-factor is determined by the sub-leading behaviour uncovered by 
perturbation theory. 

For comparison with the perturbative method, we now briefly give the exact solution 
of the model. From the rotational symmetry of the problem, only the s-wave equation is 
important. In this case the radial Schrodinger equation becomes 



The solutions of ( |2.2|) are Bessel functions, both inside and outside the well. The eigenvalue 
problem then results from the matching conditions at r = a. For a bound state, A 2 > 0, 
the matching conditions read, 




(2.1) 




(2.2) 




(2.3) 
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To draw out the similarities with superconductivity, we find the leading order solution in 
the weak coupling limit, Aa <C 1 and k <C 1, to be 

A~A exp(-^), (2.4) 

where A = (2/a) exp(-7 + 1/4). 

To make contact with the Fredholm problem, we now re-write the toy model in terms of 
an integral equation in momentum space. Using the Fourier transform of the potential, 

V(q) = ~ [ d 2 xe-^%a - Id) = -2™ 2 ^M, ( 2 . 5 ) 
a 1 J qa 

and wavefunction, 

m = J -0ji***xtpi ( 2 - 6 ) 

we get a homogeneous Fredholm equation similar to that obtained from the Schwinger- Dyson 
approach to QCD at high density, 

^--^Swf vip '- ?)m - (27) 

Partial wave analysis simplifies the solution of this integral equation. Expanding the poten- 
tial in the standard way, 

V(p-q) = J2V m (P,d)e im ^ (2.8) 

m 

V m {p,q) = 1^ d<Pe- m *V(p-q) = -^2vr f drrJ m {pr)J m {qr), 
Ztx Jo a A Jo 

we note, for a cylindrically symmetric wavefunction x(p)i that only Vq survives the integra- 
tion over the azimuthal angle in (|2.7| ) and we have 

1 r°° 

X(P) ~~ o / 9 i a9\ / d QQ v oip,Q)X{Q)- (2-9) 



2tt(p 2 + A 2 ) Jo 
We now consider the general eigenvalue problem, 

X(P) = -^A^ / W? V{f; ~ qlxiq1 ' (2 - 10) 

for which the condition A 2 = 1 identifies the specific problem of Q2.7|) . We shall find solutions, 
in the weak coupling limit, for the eigenvalues A 2 in terms of the parameters of the theory; 
k, a and A. 

Perturbatively it is convenient to work with a real and symmetric kernel. Noting that 
Vo(p, q) is already symmetric in p and q and that the kernel is real, symmetrization is 
achieved with the definition, 
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f(p) = \Jp(p 2 + A 2 )x(p), 



(2-11) 



so that the problem is reduced to the solution of the integral equation, 

f(p) = X 2 dqlC s (p,q)f(q), (2.12) 
Jo 



where 



1 / pq 

KsM = -^\J (f + A^ + A-) Vo ^ q) (2 ' 13) 

is the symmetrized kernel. 

We observe that the most divergent region of /C^ is at small p and q. Since Vq is well 
behaved at small p and q the dangerous terms arise in the denominator of ( |2.13| ) and are 
damped by A. Thus the important contributions to the integral equation arise in the region 
{p, q] ~ A. To highlight this region, we split the integral equation into two parts by 
introducing a scale 5 such that A <C 5 <C 1/a. We first look at the region p < 5, q < 5 that 
contains the leading order behaviour. In this region, Vq(p, q) ~ — k 2 7t, so the kernel may be 
split into a zeroth order term and correction: 

)C s (p,q)=K: o s (p,q)+)C' s (p,q), (2.14) 

where 



,2 



and 

.2 



**(P. = ^ (pHA ^ + A2) jo [UW) Mqr) ~ 6(6 ~ P)6{6 " q)] ■ (2 - 16) 
To the zeroth order, the integral equation ( p.!2| ) reduces to 

/°(p) = A° 2 [ S dqlC° s (p,q)f (q). (2.17) 
Jo 

We find that the only nonzero eigenvalue is given by 

_L ~^i og ^ ; (2.18) 
A° 2 2 & A' v ' 

with the corresponding (normalized) eigenfunction 



To obtain a perturbative solution to the full eigenvalue problem, we introduce a complete 
set of orthonormal functions, (jp\m) where m = 0,1,2,... labels the eigenstate and |0) 
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corresponds to the state ( [2.19|) , i.e. (p\0) = f°(p). Using the completeness relation 1 



J2m l m )( m l an d the fact that the states \m) for m > have zero eigenvalue under the zeroth 
order kernel, the eigenvalue for the full kernel can then be calculated to second order, 

y 2 = ^ + <0|X£|0) + A o2 <0|*£|0) - A o2 (0|/C^|0) 2 , (2.20) 

where the inner product is defined by, for example, 

(0\JC g \0) = j Q dqj Q dpf°( P ))C' s ( P ,q)r(q). (2.21) 

To determine the sub-leading behaviour of ( |2.20|) relative to ( |2. 18|) we are looking for 
contributions to the eigenvalue of 0(1). Since we are in the weak coupling limit, expansion 
of the Bessel functions in (|2.16|) quickly leads to the result, 

™^°fclWA))' (2 ' 22) 

so that the first order correction is found to be insignificant. This also implies that the last 
term of ( |2.20| ) can be discarded. We are then left with one term to calculate at second order. 
Although the zeroth order wavefunctions restrict the initial and final momenta to lie in the 
range [0,5], the intermediary momentum in (0|/C^|0) (resulting from the convolution of the 
kernel) is free to take all possible values [0, oo). Although no small argument expansion of 
Bessel functions with the intermediary momentum as an argument is allowed, the resulting 
integrals can be performed exactly and give, 

A o2 (0|^ 2 |0) = ^(logQ_ 7+ lV ( 2 .23) 



Substituting (|2~T8| ), ( gjp and (|2T23| ) into (|2~20|) we find that the dependence upon the 



arbitrary scale 5 that appears at zeroth order is cancelled by the second order correction, as is 
required. The leading and sub-leading behaviour of the eigenvalue has thus been calculated 
to be 



r 2 



1 k /, / 2 \ I 



A- 2 vMaJ-^iI- (2 - 24) 



From the condition A 2 = 1, the result (12.41) is recovered. 



III. THE FREDHOLM EQUATION FOR THE TRANSITION TEMPERATURE. 

We now return to the problem at hand and consider an SU (N) c color gauge field coupling 
to Nf flavors of massless quarks with the Lagrangian density, 

C = V - ^(9, - igAJty, (3.1) 
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(-PriVn) 



(-P ,-iVn) 



FIG. 1. Proper vertex function, T }' J 



c, ,c„;ci,C2 



n' ,p'\n,p). 



where F^ v = d^A a v - d u A^ + gf^A^A^ and A^ = A^t a with t a the SU(N) C generator in 
its fundamental representation. Since the Lagrangian ( [3.1|) is diagonal with respect to both 
flavor and chirality, the corresponding indices will be dropped below. 

As in [[nj, we shall study the proper vertex function, T, corresponding to the scattering 
of two quarks with zero total energy and momentum. Our conventions are the same as in BT 



but are repeated here for clarity of presentation. The Matsubara energies iv n of individual 
ingoing and outgoing quarks are labelled by n and n', respectively. Similarly, p and p' label 



the ingoing and outgoing momenta. Each of the superscripts q 



1,2, denote the color 



associated with each leg and the subscripts s, which label the states above or below the 
Dirac either + or — . Primed variables are outgoing and unprimed incoming. The 

ingoing quark pair are labelled by s = (si,s 2 ) and c = (c 1; c 2 ) and the outgoing pair by 
s' = (s[, s' 2 ) and d = (c[, c' 2 ). This labelling is illustrated in Fig. |I[ 

The proper four-fermion vertex function satisfies a Schwinger-Dyson equation (shown in 
Fig. H) which, with all indices suppressed, may be written, 



T(n' 



,P \n,p. 



T(n' 



,p \n,p 



y m 



d 3 q 

W7- 



K(n', p'\m, q)T(m, q\n,p), 



(3.2) 



where T represents the two quark irreducible vertex. The kernel has the explicit form, 

./_/__ __/ ./_ 

' s 2 (n',p'\m,q)S Sl {m\q)S S2 (-m\ - q), (3.3) 



KJ 



(n',p'\m, q 



Tf 



where S s (n\p) denotes the full quark propagator with momentum p and Matsubara energy 
i v n = ^jr( n + \)- I n order to facilitate the partial wave analysis we found it convenient 
to associate the Dirac spinors u{p) and v (p), which satisfy the Dirac equations (74) — ■ 
p)u(p) = and (74P — ■ p)v(p) = 0, to the quark-gluon vertex instead of to the quark 
propagator. Thus the vertices written in ( |3.2|) are of the form, 



Ts'-c 



U^s'^p^Usis'z, -p") r 7<5, Q /3f / a(-5l,P)^ / 3(s 2 , ~P, 



(3.4) 



with the vertex function r^^^ given by conventional Feynman rules and U(+,p) = u(p) 
or U(—,p) = v(—p), respectively. However, since the quarks are massless, 75M = — u and 
75f = — v, to simplify notation we may identify U(s,p) = u(sp). 

We complete this section of conventions and definitions with the fermion and gluon 
propagators. The zeroth order quark propagator reads, 
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FIG. 2. The Schwinger-Dyson equation. T is represented by single hashed vertices and T is 
represented by double hashed vertices. The full quark propagator is represented by a solid line and 
the bare quark propagator by a thin line. Gluon propagators, including hard thermal loops, are 
represented by curly lines. The first two orders in the expansion of T and the full quark propagator 
are given. 

S s (n\p) = % ——. (3.5) 

%v n — sp + \i 

The radiative corrections to this propagator contribute to the pre-factor C\ in ( |1.2| ). Such 
radiative corrections are depicted in Fig. |2|. This contribution has been calculated in a 
previous publication, |3l| , and the reader is directed there for details. 

We discuss the gauge invariance of this approach to the transition temperature in section 
VI, but until then we shall work in the Coulomb gauge where the gluon propagator takes 
the form, 

D u (k,oj) = D E (k,oj), (3.6) 

D ij (k,uj)=D M (k,uj)(5 ij -^), (3.7) 

k l 

D Aj {k,uj) = D j4 (k,uj) = 0. (3.8) 

In the presence of a Fermi sea, hard dense loops should be included in the gluon propagator 
in the leading order. While it is possible that a magnetic mass of order T may exist, hard 



dense loops lead to Landau damping which prevails at ^ > k^T [[35 j. In general, the 
screened electric and magnetic propagators take the form, 

DE ^)= r 2 ~ l E( r ~ V ( 3 - 9 ) 
k 2 + a E (k,uj) 

D M (k,uj) = ^ — ^ , (3.10) 

k 2 + u 2 + a M {k,uj) 

In the region k -C /i and u -C fi, both a E (k,u) and a M (k, uj) depend only upon the ratio 
x = u/\k\, x G (— oo, oo). Scaling out the Debye mass, 

by writing cr E (k, uj) = rnjjf^x) and a M (k, uj) = m|,/ M (:r), we find the dimensionless screen- 
ing functions, 
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f E (x) 



1 — xtan 1 ( — 



/m(z) = ~ 



x 

I + x 2 ) tan 1 (- ) - ./■ 

x 



(3.12) 
(3.13) 



Analysis of the functions Je{x) and fM{ x ) shows that the dangerous region for infrared 
divergences to originate in the gluon propagator is at small x, corresponding to energy and 
momentum transfer u <C k <C //. In this region the gluon self-energy functions are well 
approximated by 

f E (x) ~ 1, (3.14) 
fAf(x) ~ -\x\, (3.15) 

indicating that, while the Coulomb interaction is strongly screened by the Debye length, the 
magnetic interaction is poorly screened. It is for this reason that the attractive interaction 
responsible for the pairing instability is dominated by single collinear magnetic gluon ex- 
change. Furthermore, /e(o;) is a decreasing function of x and Jm(.x) is an increasing function 
of x, which obey the following inequalities for all x: 

f E (x) < 1, (3.16) 

f M (x) < j\x\. (3.17) 

For j > 1 we have, 

Mz) - 3^2, (3-18) 

f M (x) * |. (3.19) 

Since T corresponds to di-quark scattering it can be decomposed into irreducible repre- 
sentations of SU(N) C by either symmetrization or antisymmetrization among the initial and 
final color indices. We decompose T into its symmetric and anti-symmetric components, 

r c £(n',p'\n,p) = V2 5 c ^5^r s s , tS (n',p', \n,p) (3.20) 
+ V2 8 c ^5 c ^Tf, tS {n',p', \n,p), 

where [■ • ■] and (■ • •) denote antisymmetrization and symmetrization with weight one respec- 
tively (and likewise for T). Since the Fermi surface has a pairing instability in the presence 
of even an arbitrarily weak attractive interaction we need only focus on the attractive anti- 
symmetric channel. 

Proceeding with the partial wave analysis, we expand Tf, s (n',p'\n,p) in terms of Leg- 
endre polynomials, 

rj, tS (n',p'\n,p) = ^ l>s (n',p'\n,p)Pj(co S e). (3.21) 
j 

Using a similar expression for Tf, s (n',p'\n,p), we find from ( |3.2| ) the Fredholm equation 
satisfied by 7^ s (n' ,p'\n,p): 
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[ too 

^, (n',p'\n,p) = 7 S , (n',p'\n,p) + - ^ / dqK s , „(n',p'\m, q)^,, {m, q\n,p), (3.22) 



m,s 

where the kernel K J S , s has the form, 



K^ s (n',p'\n,p) = ^^0^ S Sl (n\p)S S2 (-n\p). (3.23) 

This integral equation can be cast into a general eigenvalue problem in the way outlined 
for the toy model in section II. To do so, we first note that the formal solution of (|3.22|) 
involves the Fredholm determinant 

V = det(l-K) = H(l-\j 2 ), (3.24) 

j 

where the eigenvalues, {A^}, are defined by the solutions of, 

fs(n,p) = - T Y, dp'K^ s (n',p'\n,p)f s ,(n',p'). (3.25) 
P 5v Jo 



Since T> appears in the denominator in the inversion of the Fredholm equation, it in fact 
governs the pole in the solution that represents the formation of a bound state of quarks 
at the transition point to the superphase. For weak coupling, we may identify the pairing 
temperature and the phase coherence temperature. 

At sufficiently high temperature all of the set of eigenvalues {A^} are greater than one, so 
that P/:0 and there is no instability — the theory is in the normal phase. As the temperature 
is reduced, we find the transition temperature to the superphase is that at which the smallest 
of {A|} reaches one. If we label the solutions so that the smallest eigenvalue is Aq, then the 
transition temperature is obtained from the inversion of the condition, Aq(Tc, g, fi) = 1. 

In order to proceed with a second order perturbative calculation, it is easiest to work with 
a symmetric kernel. In the general eigenvalue problem defined by ( |3.25| ) we have a kernel 
K s r s = K s i s > >s 8 , where the leading order behaviour only occurs for particle states (that can 
see the Fermi sea), namely for s = (+, +) and s' = (+, +). In general, however, the solution 

/++ is necessarily coupled to solutions /+_, /_+ and / by the kernels 

and K >++ , respectively. On the other hand, since ' — ' represents a quark below the Dirac 

sea, its corresponding propagator has no pole for momentum at the Fermi sea, p — fi, as 
may be seen from ( |3.5| ). In fact, the absence of this pole is sufficient to remove two orders 
of the logarithm and hence the cross-channels make no contribution to either the leading or 
the subleading behaviour. In the following we shall focus on the kernel K = K ++t++ and 
drop all particle/antiparticle notation. 

The kernel K is symmetric under simultaneous switch of the sign of the Matsubara labels 
n! and n, i.e., K(n',p'\n,p) = K(—n',p'\ —n,p). As a result, the eigenfunctions f(n,p) must 
then be either even or odd with respect to n. Since only the even functions pick up the 
singularities on the Fermi surface we shall concentrate on them and replace the sum X)^L-oo 
byE^o- 

Symmetrization of the problem proceeds by identifying 
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f(n,p) = v ^ (3.26) 



Then from ( |3.25|) and (|3.23| ) and using h{y n ,p) = h(—u n ,p) we find, 

\ 2 00 roc 

h(vn,p) = -^y. dp'K J s (ri,p'\n,p)h(v nl ,p'), (3.27) 
where the symmetrized kernel, Kg, is written, 

^(n'.j/Kp) = - J P ' - s { ^ J (n',P'\n,P)+r J (-n>,p>\n,p) 1 

This will be the starting point for the perturbative calculation of the next section. 

IV. PERTURBATIVE CALCULATION OF THE TRANSITION TEMPERATURE. 

In this section we shall concentrate upon the calculation of the transition temperature 
for a longitudinal pair in the s-wave channel, J = 0. The discussion of the transition 
temperatures corresponding to higher angular momenta longitudinal and transverse pairings 
will be deferred until section V. To zeroth order in the perturbative expansion, we start with 
the leading order solution obtained previously in f28H33|. In the following subsections we 
shall define and evaluate a perturbation series to uncover the subleading behaviour. 

For one gluon exchange, incorporating hard dense loops in the bare quark propagators, 
we find (in the K ++)++ channel), 



^=0 



n>,p'\n,p) = (l + ^) j\cos6) (4.1) 



x 



3 _ COS 0_(1 + Cos ^)|z|a_ l + costf 



+ 



{v n -v n i) 2 + \p-p'\ 2 + m 2 D f M {x) \p-p'\ 2 + m 2 D f E (x 



where x = [y n — v n i)/\p — p'\. The first term in the integrand arises from magnetic gluon 
exchange and the second from electric gluon exchange. Further discussion of this vertex 
is contained in appendix A. To obtain the leading double logarithmic behaviour we see 
that the infrared-sensitive region is in the forward direction, 9 w 0, so that cos# may be 
approximated by one in the numerator of the integrand. Furthermore, we discard the energy 
transfer term, [y n — zv) 2 , in the denominator of the magnetic propagator and approximate 
the self-energy functions, /e(x) and Jm{x), by their asymptotic forms ( |3.14j ) and ( |3.15| ) for 



small x. These approximations, prevalent in the literature, over-estimate both the electric 
and magnetic screening and thus under-estimate the attraction between quarks. How this 
correlates with the transition temperature is hard a priori to quantify. Although a weaker 
attraction should lead to a lower transition temperature and smaller condensates we shall 
see that, in fact, none of the corrections to these leading order approximations discard any 
next-to-leading behaviour of the eigenvalue. 
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With these approximations, the integration over cos# becomes trivial, 



7 J °(n,p\n' ,p') ~ 



12pp' V N 



log 



8/i 3 



4/x 2 



\p - p'\ 3 + \m 2 D \u n - v n > | + 2 l0g mf> 



(4.2) 



To leading order, we may also neglect the momentum exchange under the first logarithm. 
Although this produces a logarithmic singularity for n = n' which is potentially dangerous, 
we shall see that in fact it is well behaved. Combining the electric and magnetic parts we 
find, 



7 J °(n' ,p'\n,p) 



9 



1 



1 



where 



, H I log ■ 

I2pp'\ NJ & |z>„-zv| 



g 5 fN f y v 
256tt 4 U J 



(4.3) 



(4.4) 



Using the Ansatz, 



h°(v n ,p) 



"I + (P ~ V) 2 

and substituting into ( 3.27 ), we find that the leading order eigenvalue problem becomes 



(4.5) 



1 \ 1 x 



X^n) = Af ^ ( 1 + ^ ) - R £ i dp>K°(n>,p'\n,p), 



vl + {p- P) 2 
\ << + ip' ~ p) 



:Xj(<) ( 4 - 6 ) 



where the kernel focuses upon the most divergent region of the Matsubara energies, and we 
impose a cut, v < 5, analogous to that used for the toy model in section II: 



K°(n', p'\n,p) 



^(v 2 + (p-ri 2 )(" 2 > + (P'-ri 2 ) 



6(5 - K)0{5 - v n , 



(4.7) 



for n ^ n' and 

K°(n,p'\n, p) = 3 



to \p-p \ 



9{8-P n ). 



To calculate the leading order behaviour we move to the continuum limit, 

1 



(4.8) 



where e = 2nkBT, and perform the momentum integral for p', 



(4.9) 



1 



(4.10) 
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(assuming v n < <C Furthermore, the kernel K°, containing both — log \ v— v'\ and — log \ v+ 
u'\, can be treated using the same approximation as employed by Son [^|. This amounts 
to replacing the sum of these terms, 



log- 



1 



\v — V 



log 



1 



v + v' 



2 log 



1 



(4.11) 



where i> > = max(z>, 0'). Although this approximation may appear drastic, in fact it captures 
the behaviour remarkably well. We will show in the next section that the correction to the 
full kernel due to this approximation does not contribute a linear logarithm but that the 
continuous-discrete energy approximation does. 

In the continuum limit we incorporate ( |4.9[ ), fl4.1U| ) and ( |4.11| ) into ( |4.6[ ) , and are left 
with the integral equation, 



,2 9 



3 24tt 2 



1 \ r* du' 1 . ,, 



(4.12) 



where the ultraviolet energy cut-off 5 ~ 1 is analogous to the cut made in momentum for 
the toy model in ( |2.15| ). Note that the rescaling of v in ( |4.4| ) results in a re-scaled infrared 
cutoff, e = {Nf/2f /2 f 2 l^ . While ( [4.1 1 ) may be solved by Son's method for arbitrary cutoff 

5, the results take the simplest form when 8—1. In this case, we find the solutions 



Xi{y) = 0(5-v) 



log2/e 



sin 



1\ logl/z> 



with eigenvalues 



1 



(2j + 1) 7T 2 



'JV+1, 2 
-6AT l0g I 



log2/e 



(for 5=1). 



(4.13) 



(4.14) 



The functions hj(v,p) form a orthonormal set with respect to the continuum approximation, 
i.e. 



s du 
I 2ttJo 



dphi(v,p)hj(v,p) = Sij, 



(4.15) 



where we have dropped the contribution from the lower limit of p integration. This approx- 
imation will be adapted in the subsequent work. 

For 5^1, the eigenvalues are instead obtained by solution of the transcendental equation 
(|6.2|), resulting in 



1 

Aj 



9 



(2j + l)7T 2 



'JV+1 A 2 
-6ATl l0g I + a + 



where 



a 



-oKi + lH 



,log 3 l/<5 
log 2 2/e' 



(4.16) 



(4.17) 
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As may be seen from ( f4.9| ), the fermion Matsubara sum in the finite temperature formalism 
automatically provides a natural infrared cutoff to the continuum eigenvalue problem. Also, 
the dependence on the ultraviolet cutoff, 5, is beyond the order to which we are interested in, 
unlike in the toy model. This allows a to be ignored, so that ( |4.13|) and ( }4.14|) are sufficient 



in the subsequent calculations. More discussion on the dependence of the results on the 
ultraviolet cutoff is given in section VI. 

One concern when defining the general eigenvalue problem is that, with a non-hermitian 
kernel ( |3.23|) , in general the smallest eigenvalue defined by (|3.25| ) may be complex. Since 



the leading behavior in coupling of the transition temperature is solely determined by one 
gluon exchange, in analyzing this question we shall focus our attention on the corresponding 
eigenvalue problem with a kernel defined by ( p.23|) and (|4.1| ). In this case, by direct calcula- 



tion we have shown that the leading order of the eigenvalue is in fact real. We shall briefly 
show, further, that as long as it is not degenerate then the smallest eigenvalue must be real 
to all orders. To do so, we examine the characteristic polynomial of the eigenvalue problem, 
V(p) = detP(p) where 

To investigate the coefficients of V(p) we first consider p to be real. In this case, we find the 
relation 

Ps' A , Sl s 2 ( P ) = H^^P^sMl ( 4 - 19 ) 

by explicit examination of (|3.23| ) for the one gluon exchange in different channels, tabulated 



in appendix A together with the property that 

S Sl (n\p)S S2 (-n\p) = S* S2 (n\p)S* Sl (-n\p) (4.20) 

for a bare quark propagator. Therefore it follows that detP(p) = detP*(p), so that the 
coefficients of the characteristic polynomial are all real. Continuing p to be complex, the 
real coefficients of the characteristic polynomial imply that the eigenvalues determined by 
the characteristic equation, V(p) = 0, are either real or appear in complex conjugate pairs. 
Therefore, if the smallest eigenvalue is not degenerate it is real. Although this argument is 
specifically for the problem defined by one gluon exchange, we expect similar reasoning to 
extend to the general eigenvalue problem. 



A. The Perturbation Series. 



Given the above leading order solution, we now turn to the issue of perturbative correc- 
tions. To simplify the notation, we remove non-essential factors by rescaling the kernel and 
eigenvalues by 

K J s {n' ,p'\n,p) = (l + 1) K J s {n',p'\n,p), (4.21) 

and 
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A 2 24vr 2 V NJ A 2 

" 3 

Then to the leading order, we have Aq = Aq with 



Q 2 ( 1 , ■ 

y d , i (4_ 22 ) 



T^log?. (4.23) 



A° (J + |)7r e 

One concern, that eigenvalues of the general problem may be complex, was discussed 
in the previous sub-section where it was found that the smallest eigenvalue — the eigenvalue 
of interest in defining the transition temperature — is real if it is not degenerate. Another 
concern in defining a perturbation series with discrete energy is that the leading order 
solutions hj(u n ,p) [obtained from the substitution v — > v n in (|4.13|) 1 are not orthogonal and 
cannot be directly used as a basis for the perturbative expansion. However, each of the 
solutions may be rotated slightly so that the resulting set forms a complete orthonormal 
basis. In a more convenient notation, in which (v n ,p\j c ) — h°(u n ,p), while (u n ,p\j) = 
(jc\jc) hj(u n ,p) we can use Gram-Schmidt decomposition to write down an orthonormal 
basis, 

|J> = l/>-g=ol*>fij> , (4.24) 

where all \ j) are normalized with respect to the discrete sum over Matsubara energy, |0) = |0) 
and other \j) for j > are found inductively. The norm (j c \jc) is examined in appendix B 
along with the inner product (j c \j' c ) an d some matrix elements of the perturbative kernel. 
In this basis, we define the zeroth order operator with discrete energy, 

JCo = J2\°- 2 \i)(i\. (4.25) 

i 

This operator has exactly the eigenvalues ( |4.14| ) by definition: JCq is diagonalized by \j) and 
so \j) = \° 2 lC \j). Thus /C captures the desired leading order behavior of the full kernel at 
hand. In appendix C it is shown that in the continuum limit this operator in fact reduces 
exactly to the leading order kernel (O) that contains the double logarithmic behaviour. 



We now write the full kernel of ( ^.27[ ) in terms of this zeroth order operator and two 
correction terms, 

K J S =° = £o + Ki + /C 2 , (4.26) 

where 

K x = K° - /C (4.27) 

contains both the effects of discrete to continuous energy corrections as well as corrections 
due to the approximation ( |4.11| ) and 



/C 2 = K J S =° - K° (4.28) 
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accounts for corrections between K° and the full kernel K$ =0 . In particular, the sub- 
leading behaviour discarded by working in the continuum limit to derive ( f4.13| ) and ( [4.14| ) 
is contained in JC±. 

Following standard perturbation theory, the eigenvalue problem can be calculated up 
to second order, which proves sufficient to obtain all the next-to-leading order behaviour 
required. To define our perturbative expansion, we note that, to form a complete set, there 
are actually two types of wavefunction at zeroth order. The first type are determined by the 
kernel K° (and hence appear in the definition of /Co), are denoted by \l) (which includes |0)) 
and have eigenvalues A^~ 2 . These wavefunctions are constrained to have support only in 
the domain of K° (i.e. v e [e/2,5]), and are zero outside. The second type of wavefunction 
are not yet determined; to distinguish them from the first type they are denoted by \m) and 
have eigenvalues A~ 2 = 0. Together, these solutions give us the completeness relation, 

i = |o><o| + £lMl+£MH- ( 4 -29) 

Using second order perturbation theory and the completeness relation, we find 
1 1 



+ (0|(/C 1 + /C 2 )|0> 

o /v o 



An An 



■2 



^ (0|(/C 1 + /C 2 )|l)(Z|(/C 1 + /C 2 )|0) ^ (0|(/C 1 + /C 2 )|m)(m|(/C 1 + /C 2 )|0) 
+ 2^ i _ i +2^ 1 - n 

\2ln\ To 2 /nl/r- , ir \ln\2 



+ (0\(JCx + /C 2 )|0) + A° 2 (0|(/d + /C 2 ) 2 |0) - A^(0|(/d + /C 2 )|0) 2 (4.30) 
Aq 




__ X °o I (0|(/Ci + /C 2 )|/)(/|(/Ci + /C 2 )|0), 



which may be compared to the analogous expression for the toy model, ( |2.20|) . 

We now proceed to evaluate ( f4.30| ) and to resolve the eigenvalue problem to next-to- 



leading order in the logarithmic behaviour. This will not only determine the exponential 
and the pre-exponential factor of the transition temperature to the superphase of QCD to 
leading order in the weak coupling regime but also allow us to identify the physical origins 
of the contributions. Since the leading order is double logarithmic, the sub-leading constant 
pertaining to the logarithm is solely determined by linearly logarithmic terms. Contributions 
of 0(1) may be dropped. 



B. First order perturbation. 

The first order term in ( |4.30| ) is evaluated in two parts, the first of which is the expectation 
of K\. This itself can be separated into two parts: Kx containing corrections due to the 
approximation of the log terms given in ( |4.11| ), and JCf containing significant continuous- 
discrete energy corrections, 
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/c 1 = /cS a) + 4 b) , 

)C[ a) = 6(5 - v n )9(5 - v n >) 



(4.31] 



log 



\u —v 



+ log ,. i- - 2 log ^ 



#nn') 



3 log ft^r - 2 log ^ 



/Cf } = 9(5 - v n )6(5 



: 5 n 



/C . 



In appendix D it is shown that the contribution from /Q may be ignored 

(o|4 a) |6) = 0(1). 



(4.32) 
(4.33) 



(4.34) 



Significant continuous-discrete energy corrections do, however, arise from the expectation 
of/C? . This expression is evaluated using Zeta function techniques, and details of the 
calculation are given in appendix B. In total, we find the contribution from the first term 
to be, 



(6|/Ci|o> = -r 



7T" 



2( 7 + log 2) log- + 0(1) 
e 



(4.35) 



The second term to be evaluated at first order, K.%, contains the correction between the 
full kernel for one gluon exchange and the approximate kernel containing K° . The first order 
perturbation of ZC 2 is given by, 



(0|/C 2 |6> 



p 



2 £ 



dp / dp 



, Kg °(n',p'\n,p) — K°(n' ,p'\n,p) 



^< + {p-v) 2 )K> + {p'-v) 2 ) 



XoK)XoK') (4.36) 



s dv r s dv' 



dp / dp' 



v /( Z y2 + ( p _ yU )2)( z/ /2 + ( p /_ /i )2)- 



To produce a linear logarithmic term, or in other words to knock out one logarithmic power 
relative to the leading order, both v and v' should be ~ ksT. If either v or v' becomes large 
in comparison with ksT it will bring the integrand away from the corresponding Fermi sea 
and therefore away from the sensitive region of the gluon propagators so two powers of a 
logarithm will be eliminated. For this reason, the integrations over p and p' are dominated 
at the poles p — fi = ±iv and p' — fi = ±zz/' and the energy transfer square in the denominator 
of the magnetic gluon propagator can be disregarded. 

If the behaviour of ( [4.36[ ) is in fact linearly logarithmic it would be tied to the constant 
pertaining to the logarithms of K°. Thus ( |4.36|) demands a closer look and we start with 
the partial wave integrals of single electric and magnetic gluon interaction that combine to 
give 7 J=0 , 



2 



1 + 



N 



(4.37) 
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In the case p = p' = /x, 



7 W ~ I d C os 5 3 ~ COSd = 2 C da ^ (4 38) 

and 

«W) ~ I d cos 1 + C05g = 2 f da ^ (4 39) 

where a) = uj j ' \i and we have changed the integration variable from cos 8 to £ = sin|. We 
first note that the £ 2 terms in the numerator are free from the collinear singularity. For 
these parts the self-energy functions in the denominators can be dropped and we obtain by 
trivial integration a constant from and a constant — from ■ Amazingly, 

these sub-leading contributions of the electric and magnetic parts exactly cancel in Q4.36|) . 
As we shall in section V, this is not the case with J ^ 0. 

We now focus on the more difficult infrared sensitive parts of the above integrals, namely 

and, 

i £ 

o (441) 

For the first integral, we introduce £o such that Q <C £ *C (r/uj)^ with rj = ^3 and divide 
the integration domain < £ < 1 into the two regions; < £ < £o an d £q < £ < 1. 
Correspondingly, we have Im = I'm + Im w hh, 

^I' Vt^r (442) 



f » = f S,y + (443) 

The Landau damping approximation can then be applied to a M in (|4.42|) and the integration 
gives rise to the leading order logarithm included in K° without any additional constant. 
As for ( [4.43D , on the other hand, the mean- value theorem implies that, 

for some £ in the region < £ < £o- Since Jm{x) is a monotonically increasing function of 
x, we find, 
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Im < Co 2 /,^ < 2 f, 5 ; * 4 « 1. (4-45) 

Therefore no additional constant is found from Im either. 

The integral Ie is more subtle than Im because of the dielectric behavior at high u. In 
this case we have to introduce two cutoffs, £o and £ such that u)£o <C £ ^ ^ ^ 6o <S 1, 
and write Ie = Ie + Ie + Ie ■> where, 

^=Sl^ — , 1 2 , (4.46) 



€6 



r (c) /*o 







+ m 2 D f E { 


I) 








4// 2 e 


+ m 2 D f E 










4// 2 e 2 


+ m 2 D f E 


' £ \ 



and 

The integral gives rise to the logarithm included in K° without any additional constant. 
The mean-value theorem, applied to implies, with £ < £ < £ , that, 



j(6) = _(£q-£ )£ < £o _ Cp^L ^ i 



4// 2 £ 2 + m 2 c M|) ~ m 2 D f E (£r) e 



12 
-.0 



(4.49) 



Finally, the integral I E is divergent logarithmically at the lower limit. However this is an 
artifact of setting p = p' = fi. Without this approximation, one has to consider the partial 
wave integration simultaneously with the integrations of p and p'. The net result, with 
p — p' ~ u, amounts to a cutoff £ ~ u 4 and we expect, 

4 c )~-L^logI (4.50) 
4/i z oo 

which is completely infrared safe. Thus we have demonstrated that there is no constant 
from /C 2 + K.^ in addition to the leading logarithms in the limit to — > 0. As the result 
the first order perturbation (|4.36|) leads to no additional linear log contributions. 



C. Second order perturbation. 

Now we proceed to the last three terms of ( |4.30|) , which stand for the second order 
corrections in /Ci + /C 2 . Firstly, the penultimate term of ( |4.30| ) is easy to dispose of; although 
(0|(/Ci + /C 2 )|0) 2 ~ log 2 | this is suppressed by the factor Aq ~ (log|)~ 2 , so that it would 
contribute at most at the constant level, and not at the linear log level of interest. 

Secondly we consider the final term of ( |4.30|) , containing the sum over discrete inter- 
mediate states. For this term, the arguments given above at first order also enable us the 
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disregard /C 2 when it is sandwiched between |0) and |Z). As a result, we only need to consider 
(0|/Ci|Z) which is at most linear in logarithms. Therefore the same power counting argument 
used for the penultimate term also applies here. There is no linear logarithm associated 
with the last term and it can be dropped. 

Finally, the remaining term involves the second order matrix element 



<6|(£i 



/c 2 ) 2 |o> 



where, 



I poo poo poo 

J in - 03 H XoK)XoK') / dp J dq J dp' 



(4.51) 



(4.52) 



X 



[fCi(n,p\m, q) + /C 2 (n,p|m, g)][/Ci(m, g|n',p') + /C 2 (m, g|n',p')] 
v / K+(p-^][^ + (p'-^) 2 ] 



and 



j /*0O /"CO /"CO 



/9 J ' — ' , Jo Jo Jo 

v m >6 

[7 J=0 (n, p\m, q) + 7 J=0 (n, p\ — m, q)] [ ; j J=0 (m, q\n', p') + 7 J=0 (— m, g|n', p')} 



(4.53) 



x 



M + ( P - ri 2 ] K, + (9 - rf] K + & - rf] 



1 f,^- 1%-xMxuW) r°«v r<<9 

4 J § Z/ J £ I/' 7,5 Jo 

log + log ^77] [log pT^Tj + log ^7 



v" 2 + (q- /i) 2 

Consider J; n first. The leading logarithm of the magnetic gluon propagator has been sub- 
tracted in /Ci and we may regard /Ci as bounded in magnitude. The most dangerous con- 
tribution comes from the pole term of the integrations over p, p' and q. The summations 
over n, m and n' give rise to a contribution of order at most log 2 | (taking into account the 
inverse square root of log | from the normalization constant of each \o)- The net power of 
logarithms, when multiplied by A° 2 is thus zero. As to J ou t, since the intermediate energy v" 
is kept away from the Fermi level, the logarithms in the numerator of (|4.53|) are smeared. We 
are left with only two powers of logarithms from the lower limits of the v and v' integrations 
and one inverse power of logarithm from the normalization constant of Xo{ u )- The net power 

~ 2 

of logarithms, when multiplied by A , is — 1. Therefore we do not find any contribution to 
the linear logarithmic term from second order perturbation theory. 

In fact the only non- vanishing contribution at the linear log level, ( |4.35| ), is due to the 
discrete to continuous energy correction, and is obtained at first order in the perturbative 
expansion. Putting these results together, the expansion of the inverse square of the largest 
eigenvalue associated to the Fredholm equation with free quark propagators reads 



4 = T ^ + ^(7 + log2)log? + 0(l). 
^0 A o 



TV 



(4.54) 
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The linear logarithmic term of (|4.54| ) contributes to the factor d[ of (|1.6|) . Our analytical 



arguments on the perturbation corrections are confirmed by numerical integration, as will 
be discussed in section VI. 



V. NON-ZERO ANGULAR MOMENTUM. 

Although we have focused on the s-wave channel for longitudinal pairs of quarks, the 
collinear singularity of one-gluon exchange which is responsible for the long-range nature 
of the pairing interaction makes no distinction at leading order between zero and non-zero 
angular momenta. Therefore, in any angular momentum J channel, the largest factor con- 
tributing to the critical temperature, /c^T^, is given by the non-BCS relation (|1.3|) in all 
cases. However, the constant pertaining to the logarithm may vary with angular momen- 
tum. For the J th partial wave, the increasing number of nodes of the Legendre polynomial 
Pj(cos9) with increasing J reduces the effective attraction between quarks. As we shall see, 
although formally of 0(1), the suppression for J ^ can be of several orders of magnitude. 

To explore such a correlation, we evaluate the J ^ partial wave component of the 
irreducible vertex function, 

7>VKP) = -j (l + ^) (2J+ 1) j\(cos9)Pj(coa9) (5.1) 

3-cos#-(l + cos#)|Ef^ l + cos£ 



V n - IV 



_l_ \p- p'\ 2 + m 2 D f M (x) \p-p'\ 2 + m 2 D f E (x) 



To uncover the dependence on J, we set p = p' = fi. The numerator of the integrand can be 
decomposed into a constant term and a non-constant term with the latter vanishing in the 
forward direction, cos# = 1. Since Pj(l) = 1 for any J, the contribution from the constant 
term is independent of J and can be approximated by the J = vertex given by Eqn. Q4.2|) . 
For the non-constant part of the numerator, proportional to 1 — cos#, the forward direction is 
infrared safe, and the self energies /a/ and Je may be dropped to the order we are interested 
in. The resultant integral produces a sub-leading J-dependent term. We find, 



r J (n',p'\n,p) = --|^ (l + 1) (2 J + l)(log 1 — + 3 CJ ), (5.2) 



with, 



f 1 , Q Pj(cos9)-l 

cj= dcos9— — . (5.3) 

J-i 1 — cos 9 

Using the recursion formula for Legendre polynomials, 

(J + l)P J+1 (cos#) - (2J + l)cos#Pj(cos#) + JPj_i(cos#) = 0, (5.4) 
we find that, 

(J + l)c J+1 - (2J + l)cj + Jcj. x = -25 J0 . (5.5) 
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Since Co = 0, the solution to Q5.5Q for J > 1 is, 

cj = ~2Y<l- (5-6) 

n=l ' t 

Following the arguments of the perturbation theory for J = 0, we see that cj contributes 
only to the first order perturbation of the eigenvalue, and hence the critical temperature in 
the angular momentum J channel reads, 

k B T^° = e^k B Tf (5.7) 
~ 2.5 x KT 3 k B T^ =0 for J = 1, 

where Tq =0 is the corresponding critical temperature in the J = channel. 

The foregoing consideration can be generalized to transverse pairings between opposite 
helicities. The only modifications include replacing 3 — cos 9 in ( |5.1|) by 1 +cos 9 and Pj (cos 9) 
by d(i(9), the Wigner D-functions ]3B]. The final formula for the transition temperature is 



given by (|5.7p with cj replaced by 

If J J + 1 \ . . 

2^+2jTI C ^ 1 + 2JTl Cj - 1 > (5 ' 8) 

with J > 1 and Cq = 0. Relative to the longitudinal pairing at J = 0, we find that for 
J = 1 a longitudinal pairing is suppressed by a factor of e -6 whereas a transverse pairing is 
suppressed by e~ 4,5 . 

Physically, ( |5.7| ) shows a large suppression of the transition temperature with non-zero 
angular momentum relative to that for s-wave pairing. As the quark-gluon plasma is cooled 
from the normal phase, the first critical temperature reached is Tq =0 and the corresponding 
long range order is then switched on. This would then invalidate the relation (|5.7|) for J ^ 0, 
which was derived using normal phase propagators. Furthermore, there is no simple relation 
between Tq^° and the corresponding gap energy at T = 0. Nevertheless, the suppression in 
( |5.7|) indicates that if for some physical reasons the long range order with J = cannot be 
formed, then the J ^ pairing can only happen at very low temperatures. 



VI. NUMERICAL RESULTS AND GAUGE INVARIANCE 

In the previous section we have examined the perturbative expansion of the full eigenvalue 
problem, ( |3.28| ), and have found a single correction term at the linear log level, arising from 
the discrete to continuum energy approximation, ( |4.35|) . While this result is simply stated, it 
came from an involved perturbative treatment of the Fredholm integral equation to second 
order. Thus, as a check, the analytical arguments on the perturbation corrections have 
been confirmed by numerical evaluation of both the first and second order terms in the 



perturbative expansion, ( 4.3C ) 



At this point, a comment on the integration limits is in order. We recall that a cutoff, 
v < 5, was imposed in order to control the zeroth order kernel, ( [4.7| ). As this cutoff is not 



otherwise physical, its dependence ought to drop out of the final result. This in fact has been 
demonstrated in the toy model of section II, where the cutoff dependence cancels between 
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the zeroth and first order terms. For the actual case at hand, however, the cutoff dependence 
behaves differently. Consider the zeroth order one-gluon exchange kernel, ( [4.12 ), which may 
be rewritten in terms of the variable x = log v: 



f(x) = -A 2 / dx'x > f{x') 



(6.1) 



where a = log | and b = log 8. The exact eigenvalues of (|6.1|) may be obtained using Son's 
method ||28|| . The resulting spectrum is given by the zeros of the transcendental equation, 

6A = -cotA(6-a), (6.2) 

and may be developed perturbatively for (b/a) <C 1. We find, 



7T 



log 



m 2 



Ml 



+ o 



log 8 



(6.3) 



which corresponds to Eqn. ([4.16|) for 8 ^ 1. This indicates that the cutoff dependence of the 
zeroth order term is in fact down by three powers of the logarithm, and hence contributes 
below the level we are interested in. This in fact explains why, unlike the toy model, we 
find no cutoff dependence at the linear log level term by term in the perturbative expansion 
( |4.30| ). Thus our numerical work for the above was performed with a fixed upper cutoff of 
8=1. 

In addition to numerical confirmation of the perturbative expansion, we have also inves- 
tigated the numerical solution to the Fredholm equation ( 3.28|) . This provides both a direct 
verification of our analytic results and a preliminary demonstration of the gauge invariance 
of our approach. Working in the normal phase makes treatment of gauge invariance easier 
to handle. In particular, as discussed in |3~I] , the Fredholm determinant can be represented 
as a series of bubble diagrams, the sum of which is gauge invariant. 

As this numerical work is secondary to the analytical results presented above, we make 
several simplifying approximations. Firstly, to avoid infinite sums, we only work with a 
continuous energy integral. As a consequence, we will not reproduce the discrete to con- 
tinuum energy correction of Eqn. ( f4.35| ). Secondly, we a priori dispose of the momentum 
integral by noting that the dominant contribution to the kernel occurs near the Fermi sur- 
face, {p, p'} ~ ji. These two approximations correspond to the steps taken in Eqns. ( fO| ) and 
( |4.10| ) in deriving the zeroth order integral equation, ( |4.12j ). Numerically, however, we do 
not employ Son's approximation, ( f4.11p , and instead use the complete one gluon exchange 
vertex (evaluated with momenta at the Fermi surface). 

One difficulty immediately arises as a consequence of setting p = p' = fi in the one gluon 
exchange kernel, and that is the spurious divergence in the Coulomb propagator, ( |3.9| ), in 
the forward direction. While this divergence is eliminated by a full treatment incorporating 
momenta off of the Fermi surface, as shown in ( 4.50 ), here we have no such luxury without 
reintroducing the momentum integral. This difficulty prompts us to work in the covariant 
gauge, using the propagators of p9[. Using the covariant gauge propagators directly in 



numerical analysis of the integral equation provides a first check upon the gauge invariance 
of the approach. 
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FIG. 3. The largest eigenvalue, 1/A 2 , of the Fredholm problem (|6.4j ) as a function of log ~. The 
numerical solution is given by the solid line, while the zeroth order solution, ([T7]), is given by the 
dashed line. 



Taking everything into consideration, we thus examine the numerical solution to the 
integral equation, 



s dv' 



with kernel 



(6.4) 



(6.5) 



For the one gluon exchange interaction in covariant gauge, we use the result of [J29| . Namely, 
in our notation, 



7(1/,!/) = 2/ d& 



(1/ - z/) 2 + Ve 2 + m 2 D f M (x) (y - u'f + 4/i 2 e 2 + m 2 D (l + x 2 )f E (x 

(6.6) 



where x = \v — v'\j (2/i£). Note that the electric gluon exchange term is modified in covariant 
gauge, as may be seen by comparison with ( 4.38Q and ( f4.39|) . For this linear Fredholm 
problem, which is similar to the non-linear gap equation, the numerical solution for the 
largest eigenvalue is easily obtained by iteration. One subtlety arises in this case, however, 
in that the relation between Tq and the gauge coupling g is given by inversion of the critical 
temperature condition, A 2 (T C , g, fi) — 1. To investigate the behavior of (|6.4| ) near the critical 
point, we thus use the zeroth order relation, g 2 (l + jj) = 67r 4 (log - )~ 2 , where e = 2nkBTc, to 
fix m 2 D in ( p.6[ ). For explicit numerical work, we furthermore set 5 = 10 10 fi, thus employing 
a finite but large cutoff. 
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log(2/e) 

FIG. 4. The difference between the numerical solution and the analytic result, = ^ k>g 2 f 

(solid line). For comparison, the ^-dependent behavior of the eigenvalue, given by ( |6.2| ), is indicated 
by the dashed line. Note the different scale used on the vertical axis, compared with Fig. |3[ 

The numerical solution to fl6.4|) is presented in Fig. |3], and shows good agreement with 
the analytic result, 

1 4 , 2 2 4 / 2 , 256/i 6 \ 2 
TT = — log 2 - = — log-+log ^ , 6.7 

at least for log - ^> 1, corresponding to g <C 1. This behavior is highlighted in Fig. [|, where 
the difference between the numerical and analytical solutions is shown. Up to the linear 
log level, there is no discrepancy between the solutions; as they were obtained in different 
gauges, these results indeed support the gauge invariance of the critical temperature. 

The behavior in Fig. [| clearly indicates the presence of subleading terms, contributing 
below the linear log level. A large portion of this appears to be related to the treatment of 
the ultraviolet cutoff, taken as 5 = 1 in the analytical work. The numerical work, on the 
other hand, is suggestive of the presence of an effective cutoff near the Fermi surface, S ss fi, 
corresponding to 5 = (it / 256) (m p / fi) 5 . While any difference between these cases is formally 
seen from ( |6.3|) to be down by three powers of the logarithm, it however contains a hidden 
enhancement due to the factor of m b D entering into 5. Since mu ~ g ~ (log (using the 
zeroth order relation), the relevant expansion parameter is (log 5/ log 2 ) ~ (51ogx)/x where 
x = log -. Even for x ~ 100, which is the largest value we examined numerically, its cube is 
only down by about one percent. 

For relevant values of g, which are of order one, the values of x range from about 5 to 
30, in which case this sensitivity to the cutoff can lead to effects of order 10%. On the 
other hand, for the nonlinear gap equation, this sensitivity appears much more pronounced; 
varying the ultraviolet cutoff between 5 and /i8 can lead to as much as a factor of two or 
more change in the value of the zero temperature gap! 
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VII. CONCLUSION. 



In this paper, we present an exact calculation of the transition temperature to the color- 
superconducting phase of SU(N) C QCD in the limit of high baryon density. To do this 
we cast the Fredholm integral equation that arises from the Dyson-Schwinger approach to 
di-quark scattering into a general eigenvalue problem. The leading order of the eigenvalue, 
found to be doubly logarithmic, gives rise to the non-BCS dependence of the exponent upon 
the coupling. We have developed a perturbation method which enables us to calculate the 
eigenvalue beyond the leading order in logarithm. The result, together with the contribution 
from the radiative correction of the quark propagator found in [BTJ , leads to a complete 



determination of both the exponential and the pre-factor of the scaling formula for the 
transition temperature, (|1.2j ), to the leading order of the running coupling constant. The 
transition temperature and the zero temperature gap are both suppressed relative to previous 
claims due to the effects of the quark self-energy. For N = Nf = 3, these effects reduce the 
results by a significant factor of 0.17663. 

Using the one-loop formula for the running coupling constant at the chemical potential, 

9 127T 2 , , 

9 2 = ,„„ -tt - , (7-1) 



%N-N f ) logf 

the behavior of the transition temperature is shown in Fig. |5] for a longitudinal pairing with 
zero angular momentum for N = Nf = 3 and A = 200MeV. We note that g — 1 corresponds 
to \i ~ 10 6 MeV and a transition temperature ksTc — 1.9MeV. Our calculations become 
valid in the limit of high baryon density where the theory is weakly coupled, corresponding 
to the upward sloping arm on the right of Fig. ^ where \x > 10 6 MeV. However, there is 



some contention that such perturbative series are valid for g < 4, see [|33| and the references 
therein. Extrapolating our results to g ~ 3 corresponds to fi ~ 500MeV and ksTc — 
3.7MeV. This is the point at which the graph in Fig. ^ turns over and the behaviour becomes 
unphysical. This limit suggests that for g > 3 strong coupling effects truly dominate and 
the transition temperature is determined by the mechanisms of chiral symmetry breaking 
and confinement but for g < 3 our results may prove increasingly accurate. 

The ratio between the transition temperature and the energy gap at zero temperature, 
A, is an important number for superconductivity. This ratio has been calculated for zero 

dentical to that of BCS theory, i.e. 

vre" 7 . (7.2) 



angular momentum in |HJ and the result is identical to that of BCS theory, i.e. 

A 



k B T c 

Combining our formula for ksT(l =0 , (|1.2j ), and this ratio we obtain that, 

/ O \ 2 / 6 AT 7T 2 



A = 2c;A = 512cX I — J ^e"V -+ 1 • , (7.3) 

with c[ given by ( |1.5| ) and A given by ( |1.3| ) and previously determined in |29[. Apart 
from the new factor c[, this is in agreement with the gap energy presented in |33[ (we have 



highlighted the factor of 2 which was found numerically in |29[ and by analytic means in 
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LI (MeV) 

FIG. 5. The transition temperature, ksTc, as a function of the chemical potential, fx, for 
N = Nf = 3 and A = 200MeV. 

||33||). What is significant, and perhaps somewhat surprising, is that the additional unknown 
pre-factor of 0(1) present in all previous calculations of the zero temperature gap has been 
found to be exactly one to the leading order in the coupling. 

We have extended the analysis to find the transition temperature for both longitudinal 
and transverse pairings for arbitrary angular momentum, J > 1. Relative to the longitu- 
dinal J = pairing these channels are suppressed by several orders of magnitude. As the 
temperature is reduced, once the highest transition temperature is reached the results for 
lower transition temperatures are invalid. However, if for some physical reason the J = 
longitudinal pair cannot form these results indicate that color-superconductivity may only 
occur at very low temperatures. Further investigation of the ratio (|7.2|) generalised to J > 
may produce interesting results for condensates of non-zero angular momentum. 
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APPENDIX A: SINGLE GLUON EXCHANGE VERTEX. 

In this appendix, we set out some conventions and tabulate the single-gluon exchange 
vertex in all channels with s = ± and s' — ±. In Coulomb gauge, the relevant one gluon 
exchange contribution to ( |3.21| ) reads 



g 2 1 
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+ -jD^(A;,a;)[ W t( s ' 1 ^)«(«iP)]K(-4p / )«(-^)], 



(Al) 



where k = p' — p and uj = v' — v. Using the Dirac equation (— ij ■ p + 74p)u(p) = and 
introducing the notation 



and 



^s'l! s 2' s l> s 2 



[u{s' l p')^u(s 1 p)][u(-s' 2 p')^u(-s 2 p)} 



(A2) 



we have 



S ' s 1 ,s 2 ;s 1 ,s 2 



u^(s' 1 p')u(s 1 p)][u t (-s' 2 p')u(-s 2 p)} 



(A3) 



Tf, s (n>,p>\n,p) = ( 1 + - ) D M (k,u 



(s'-j/ - stfXs'zjJ - s 2 p) 

V s' 1 , S ' 2 ;s 1 ,s 2 ,^ t _^ 2 O s > i>s > 2 . s ^ S2 



l + -)D*(k,u)S s[>s , 2 ., Sv 



(A4) 



Our convention for gamma matrices is that {7^,7^} = 5^ with each 7^ hermitian and 



75 = 71727374 



-I 
-I 



(A5) 



The coordinate system is chosen such that p\\z and p' _L y with p ' ■ x > 0. The azimuthal 
angle of p is assigned to be and that of — p to be ir. The polar angles of ±p and ±p' are 
p : (0, 0), — p : (ir, 7r), p' : (6 1 , 0) and — p' : (n — 9,ir). The chiral spinors at ±p and ±p' are 
related via 



u(q) = e _ i CT3 V 



2 cr 2 V e2 (T3' 



(A6) 



with q any one of — p and ±p' with (#, </>) the corresponding polar angles. For 75«(p) = —u(p) 
we write 




1 



(A7) 



Then the expressions for V s > s > - s s and S s > s > - s s follow and are listed in Table. I. 



APPENDIX B: FIRST ORDER CORRECTIONS FROM /cS B) . 

In this appendix, we shall calculate the normalization constant of \j) due to discrete- 
continuous energy corrections, and calculate the contribution of Kx giving rise to ( [4.35 ). 
We begin with a simplification of the matrix elements of an operator of the structure 

y, P '\o\v,P) = Mv ') ] (bi) 

v 2 + (p — iifdv' 2 + (p' — jj) 2 
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s l s> 2i 


s 1 s 2 




^ s l i s 2 ' s l ' s 2 


+ + . 


+ + 


3 
2 


— 2 COS (7 


5 + 5 cos 


+ -, 


+ + 




2 bill (7 


— ^ sin 9 


- + , 


+ + 




1 cin ft 
— 2 ^m C 


2 1 bin# 





+ + 


3 
2 


1 m« # 
— 2 CUo (7 


— 1 + 5 cos# 


++, 


+ - 




— i sin 9 


2 1 sin# 


+ -, 


+ - 


1 

2 


— 2 COb (7 


5 + 5 cos 


- + , 


+ - 


1 

2 


4_ 1 cc\c ft 
t 2^ cub u 


\ — \ cos 





+ - 




2 bin u 


2 1 sin# 


++, 


- + 




1 •„ /i 

2 Sin t7 


— \ sin 


+ -, 


- + 


1 

2 


+ | cos 9 


i — 1 cos 


- + , 


- + 


1 

2 


+ \ cos 9 


5 + 5 cos 




- + 




— \ sin 


1 • n 

— 5 sintf 


+ + , 




3 
2 


— \ cos # 


— I + 5 cos# 


+ -, 






2 1 sin 9 


— \ sin 


- + , 






— \ sin 9 


\ s'm.9 






3 
2 


— \ cos # 


5 + 5 cos 



TABLE I. The factors V and S for all combinations of incoming and outgoing quarks above 
and below the Fermi Sea. 

For a function (v,p\ip) pertaining to the ket vector we introduce the function 

J ° yV 2 + (p - fx) 2 



so that 



(ip'\0\ip) = jr 7 - Y7 -y(^> n /)o(z/ n /,z/ n )(z/ n |^). (B3) 

„=o (n' + i) (n + |) 

Note that we have extended the lower limit of the p and p 1 integrals to —00; the correction 
introduced is beyond the order of interest. The same simplification extends to the next 
appendix as well. 

Given Ki in Eqn. (|4.33 ), we wish to evaluate the contribution 



_ , tt _ (O c |log^|0 c ) 1 

^ w= w (B4) 



Note that while the eigenf unction, 



MOc) = xoiy) = 2 V^ sin ( tlo &l ) • (Br)) 



where 
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' = 2l|| « X ' (B6) 

is properly normalized in the continuum norm, this may no long be the case under the 
discrete norm. Thus we must include the denominator (0 C |0 C ) in (|B4[), which is to be 
calculated with discrete energy variables. 

Using (|B5|) , this discrete norm of the zeroth order eigenfunction reads, 

At N ° e / 1 \ 

(0 c |0 c ) = -]T^sin 2 Hog- , (B7) 

7T ^ U n V V n ) 

where 

v n = (n + |)e (B8) 

and iVo ^> 1 is the discrete cutoff corresponding to the ultraviolet cutoff, S, in the continuum: 
5 = (No + \)e. 



Evaluation of this sum may be performed in terms of the generalized zeta-function, |37 



OO 1 

C( J . fl ) = Em;' ( B9 ) 



resulting in, 



(o c |o c > = |{c(i + o + , i) - C(i + o + , N + |) 



-Rer 2tt 

Using the asymptotic formula for t — > 



C(l + 2it, |) - C(l + 2it, N Q + §)] }. (BIO) 



C(l + 2zt, a) = ^1 ... + 0(t), (Bll) 

2tl [a + it) 



which follows the integral representation of ((s, a), we obtain that 



C(l + 2it i)~-il-V>(^), (B12) 



2 7 2t \2 
with ^(z) the logarithmic derivative of T(z) and 



((l + 2tt,No + h~-UN + l) 2t \ (B13) 



2' 2t 
Therefore, we end up with 

At 2 
(0 C |0 C ) ~ 1 - -Ml) + log 2) = 1 + 2( 7 + log 2) log" 1 — , (B14) 
7r e 

where we have used = — 7 — 2 log 2. 

Evaluating the matrix element of the kernel is more difficult. We start with the expression 
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No i N 



(O c | log ^|0 C > = - £ 4- E 4- sin ( t log 4- 



71 n =0 Vn m=0 V ™ 









(t log 


log^ 


sin ft log -7—^ 




Z/ >. 





(B15) 



In order to resolve z>> we split the double sum, making use of the symmetry under n «-> m 
interchange, and obtain 



1 



it 



^0 f d / 1 \ 

(0 C | log ^|0 C ) = --E —S n (t)- cos t log — 



(B16) 



n=0 



where 



S'n(t) = E — sin ft log — ) = Im {r i4 [C(l + it, \) - C(l + it,n + §)]}. (B17) 



m=0 u m 



Inserting ( B17|) into ( |B 1 6|) , we find 

e-'C(l + 2t,|)-e l *C(l-zt,i; 



. h 1 . , lit 
(0 C | log — |0 C ) = — 

V > 7T 



(B18) 



d_ 
dt 



i~ u (ai + it, \) - C(l + it, N + §)) + e u (C(l - it, i) - C(l - it, N Q + § 



_ r it 



di 

dt 



r u (z(t',t; i) - Z(t!,t; N + §)) + e*(Z(t', -t; §) - Zit! , -t; iV + j 
e-" (Z(-f , t; |) - Z{-t', t; N Q + §)) + e u (Z(-t', -t; \) - Z{-t', -t; iV + §] 



t'=t 



where, 



Z(t / ,t;a)^E (n+ ; )Wt C(l + /VM + ( / + ;» 



(B19) 



r(i + if)r(i + it) 



/■OO /"CO 

Jo dx L dy JT^ 



(1 - e _ *)(l - e~ x -y) 



Changing the integration variables to u = 1 — e x and v — 1 — e x y , we find for t — > and 
t' -»• that 



Z(t',t;a) 



r(a) 



t'(t + t') T(a + i(t + t')) 



Therefore 



and 



Z(t',t; 



Z(t',t;iV + -) 



1 | 1 

't'(t + t') + t'^V 



. 1 (iv + 
t'(t + t'y t 



(B20) 



(B21) 



(B22) 
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Upon substituting the above expressions into 
manipulations we obtain 



and expanding for t <C 1, after some 



(0 C | log ^|0 C ) ~ i[l - -(V>(§) + log 2)] = 3 log 2 H[i + 4(7 + log 2) log" 1 -]. (B23) 



7T 



7T 



Finally, after combining this with the discrete normalization, ( |B14| ), and subtracting out 
the zeroth order eigenvalue, we find the claimed result, 



(6|/cf)|0) = - 2( 7 + log 2) log; + 0(1) 



4 



7T 



(B24) 



In addition, we have also calculated the general inner products (j c \j' c ) and the matrix 



element (0 C | log^|j c ), and found that 



(Mjc) = (0 C |0 C ), 



(j c \j' c ) = log 



for j ^ f, 



and 



(O c |log^-U c ) = o(log|) 



(B25) 
(B26) 

(B27) 



for j 7^ 0. 



APPENDIX C: THE DISCRETE KERNEL IN THE CONTINUUM LIMIT 

Here we demonstrate that in the continuum limit /Co indeed leads to the continuum 
approximate form (|4.11|) . To leading order, there is no distinction between \j), j c ) and \j), 
nor between the eigenvalues of the continuous or discrete theories. Thus we are interested 
in the quantity, 

(v„|fcoM = Y,K*(vn\3)QK>)^Y.^\v\j)W)i ( C1 ) 

3 3 

where the last expression on the right hand side contains continuous energy variables and 
the last inequality holds at leading order only. Since (v\j) = x( u )i using equation ( |4.13|) , we 
may re-write (|C1|) as, 



E A r 2 (^b')0>') = -2 lo s ? w» - M - + 6 ")\ > ( C2 ) 



where 



oo -y 

J («)=Et -2C0S[(J + |)«], -7T<a<7T, (C3) 

3=0 [j + \ 



and 
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0„ = 7T 



log 1/0 



log2/e' 

To evaluate 1(a) we consider its derivative, 



o < e u < it. 



(C4) 



rfa J + 2 



e i(i+5)a 



-Imlog 



1 + e 



jq/2 



1 _ gM*/2 



2m7r , 

2' 



(C5) 

where m is an integer, yet to be determined. The condition 4^- = — £ then fixes m to 
be zero and I (a = it) = fixes the constant of integration to give 1(a) from ( |C5| ) . Since 
1(a) is an even and antiperiodic function, we find, 



Ka) 



71 



7T 



M) 



-2tt < a < 2tt. 



(C6) 



Substituting flCTf ) into flC2|) then yields, 

{loe l/z>' if z>' > z>, 
/ 
log l/z> if z>' < 0, 

which is in direct agreement with the leading order kernel written in ( |4.11 ). 



(C7) 



APPENDIX D: FIRST ORDER CORRECTIONS FROM K,[ A) . 

In this appendix, we calculate the first order perturbation from . On splitting the 
discrete energy sums resulting from evaluation of the matrix element into diagonal and 
off-diagonal parts, 

(0|4 a) |0) = / + J, (Dl) 

we have 

/ = ^E/ d P dj/KT\n,p\n,p f )(0\n,p)(n,p f \0) (D2) 

n=0 J ° J ° 



and 



J =^EE(!-U dP dp'IC^(n,p\n',p')(b\n,p)(n',p'\0), (D3) 

P n=0n'=0 J ° J ° 



where Nq 1, the discrete cutoff, is given by 5 = (No + ~)e with S the ultraviolet cutoff 
introduced in section IV. 

For the first term, using ( 4~8|) , we find that 

St No r°° r oo 3 log f% - 2 log i , 1 N 

' E f dp I dp' j= , g - = sin 2 (t log 1 , (D4) 



Vn 
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with t = r) . 7T -2 . For e <C u, 

21og 7 ^' 



2t ~ 1 
" ^ 3 k (« + I) 2 



sin' 



can be approximated by 



(tlog— j / da:/ dy 

V V n > J-oo J-oo (X Z + 1)( 



l)(y 2 + l) 



x 



3 log 



2/i 



x - 2/|(n + |)e 



21og- 



n+|)e 



2*, 2/i ~ 1 , 1 

-—log— — -^sin (tlog — 

= 0{l) (D5) 

for t — > 0, since log ^ ~ and the infinite series converges. 

Next, instead of calculating J directly, we evaluate the more general expression, 

X N go 

J m = -as E E(* _ s nn') / dp/ dp%f\n,p\n\p')(m\n,p){n',p'\ih), (D6) 
so that J = J m =o- It follows from (|4.32|) for e <^ fi, that 



2t 



N N 



7T ±^ n (n + |)(n' + i) v i/„ - i/„ 



+ log- 



n=0n'=0 V" 1 2^'" 1 2 

x sin 

No n-1 

7rS^,(n+i)(n' + i)^ZVi> n 
4t ~ 1 dz> dz>' /z>\2/ 



21og^ 



1 1 r 1 
(2m + l)tlog — sin (2m + l)t log 



£1 jvo n-i °° 1 / 77 / \ 27 r 1 1 

— EE? — i iw rr i\ Et(^ sin (2m + 1)* log — sin (2m + 1)* log — 

ir * — ! id -J- inn.' -I- — I r — ' /. \ i/_ / L j> J L r/_/ 



7/ J | V' 



sm 



(2m + l)i log — sin (2m + l)tlog — 



7/ 



OO -I /-log - /•£ 

— E T / ' ^ / dx'e" 2 ^"^ sin[(2m + sin[(2m + l)ta'] 

7T l=1 Z Jlogi Jlogi 

Aj. OO 1 /-log 2 fX 

— E T / ' / d^V 2 ^"^ sin [(2m + l)te] sin [(2m + l)ta'l 
I I Jo Jo 



1 OO 

~-E 



7T 



2^Z 2 + (m+±) 2 t 2 12 



(D7) 



as t — > 0. Therefore J = 0(1) and no linearly logarithmic terms result from (0|/C^|0) = 
(9(1). On the other hand, as m — > oo at a fixed £, J m ~ 1/m. Since J m is proportional to 
the first order shift of the -4- given by ( |4.14| ), the large m behavior of J m explains the weak 
logarithmic singularity of the Fredholm kernel. 
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